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Abstract 

We calculate the shear viscosity to entropy density ratio in pure SU(N) Yang-Mills theory 
below the critical temperature using the strong coupling expansion. The result for the n/s 
ratio for temperatures around the phase transition is 0.22 for N = 2 and 0.028./V 2 for N > 2. 
The results are consistent with the conjectured 1/4-ir lower bound inspired by the AdS/CFT 
correspondence. 

1 Introduction 

Ideas originating from the AdS/CFT correspondence found a surprisingly rich set of applications 
in the description of strongly coupled QCD plasma. At first sight this may indeed seem surprising 
since the prototype case of Af — 4 supersymmetric Yang-Mills theory is very different from QCD. 
For instance the former is a conformal theory while the latter is not. However, the application is to 
QCD in the dcconfined phase and for large temperatures QCD is actually close to being conformal. 
There are other reasons too why it is perhaps not an accident that AdS/CFT inspired ideas work 
better than expected [1]. 

One particular instance where AdS/CFT based arguments provide a good description of experi- 
mental data is the case of transport coefficients. In heavy ion experiments a fairly good description 
of the data can be achieved using (classical) hydrodynamical models and the coefficients of the 
corresponding Navier-Stokes equation can be determined from measurements. These and in par- 
ticular the V2 value of the anisotropy of the observed flow [2, 3] point to the fact that the QCD 
plasma is an almost perfect fluid. A dimensionless quantity can be formed from another coefficient, 
the shear viscosity rj and the entropy density s which experimentally turns out to be rather small, 
r)/s w 0.2 [1,4]. 

Based on the aforementioned AdS/CFT ideas a universal lower bound was conjectured rj/s > 
1/4-7T for a large class of quantum field theories [5-7]. This bound is known to be saturated by 
Af = 4 super Yang-Mills theory at large 't Hooft coupling and large TV in which limit the gravity 
dual is weakly coupled and the calculation is tractable on the gravity side. The first corrections at 
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finite N and 't Hooft coupling was shown to be positive and hence comply with the conjectured 
bound [8-11]. 

One argument supporting the universal lower bound conjecture is that the energy-time uncer- 
tainty relation can be reformulated entirely in terms of rj/s using the quasi particle picture [6]. The 
lower bound of the uncertainty relation yields a lower bound for rj/s. 

However the precise set of quantum field theories where the bound is expected to hold is not 
clearly formulated. Counter examples are known which served to limit the conjectured validity of 
the bound [12-14]. As a result it is not clear what precise statement should be proved or disproved 
and a case-by-case approach seems more fruitful [15-18]. Some of the studied cases satisfy the 
bound while some do not. 

The counter examples in [12, 19] fall into two types. In one of them the number of particle 
species is very large, in the other the number of excitations of a particle are very large. While these 
constructions are convincing perhaps the conjectured bound can be valid in a limited sense. In any 
case this issue is far from being settled [13, 14]. 

The second type of construction in [19] may be realised in the framework of quantum field theory 
as infinitely many peaks finally melting into a continuum. This suggests that a good candidate for 
violating the proposed lower bound should contain a strong continuum preferably with no quasi 
particle peaks in the spectral function i.e. it should be neither Bose nor Fermi gas. 

Apart from models with known gravity duals analytic computations were done in QCD using the 
quasi particle approach [20-22] and recently using perturbation theory in QED [23] . These studies 
yield rj/s ~ 1/g 4 In g^ 1 . Clearly, when perturbation theory is applicable the coupling is small and 
the above ratio is large. And vica versa, a small measured value of this quantity hints at a large 
value of the coupling constant. Newer calculations [24,25] suggest however that higher order quasi 
particle processes may still lead to small viscosity even in the perturbative regime. 

A non-pcrturbative calculation of rj/s has been attempted using a lattice discretization in SU(3) 
pure gauge theory [26-28] . Unfortunately even in cases where a reasonable ansatz is available for the 
spectral function e.g. a transport peak structure, quantitative details are very difficult to extract 
from correlation functions [29] . 

Using the large Nf expansion of the 2PI effective action, which formalism sums an infinite set of 
diagrams, transport coefficients were calculated in [30] . It was found that for large fermion masses 
the shear viscosity can be arbitrarily small. 

In this paper the rj/s ratio of pure SU(N) Yang-Mills theory will be calculated on a Euclidean 
lattice at strong coupling for T < T c . The actual final result is 0.22 for N = 2 and 0.028.ZV 2 for 
TV > 2 at temperatures close to the deconfincment transition. We find a vanishing rj/s ratio at 
T — > 0, however the systematic uncertainties inherent in the strong coupling expansion are hard to 
estimate. Similar calculations for transport coefficients have been carried out in [31,32]. 

The strong coupling expansion has several conceptual difficulties not specific to the observable 
we calculate [33]. Since on the lattice with decreasing lattice spacing the bare coupling constant 
is going to zero, we cannot approach the continuum limit. The commonly used strategy is trying 
to express the results in terms of physical quantities, which are then interpreted as renormalizcd 
quantities. Another issue is the presence of the roughening transition which makes it impossible to 
reach temperatures higher than a certain value. Despite of these difficulties, however, the strong 
coupling expansion not only works rather well for describing several concepts in the confined phase 
such as chiral symmetry breaking and confinement but also quantitatively [33] . 

The outline of the paper is as follows. First the shear viscosity is calculated in section 2. The 
calculation is based on the Kubo formula, which will be related to the lattice plaquette-plaquette 
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correlator. In the strong coupling regime the plaquette-plaquette correlators are dominated by 
elongated tube-like surfaces. A subclass of such simple surfaces as well as more complicated ones 
involving a loop wrapping around the finite imaginary time direction are resummed in section 2.6 
leading to a final formula for r\. In section 3 the entropy is calculated from the derivative of the free 
energy. In section 4 we calculate the shear viscosity to entropy ratio in the limit of the temperature 
being close to the transition point. Section 5 concludes with a discussion of our findings. 



2 Shear viscosity 

The shear viscosity is one of the linear response coefficients coming from the response of the energy 
momentum tensor to an external current coupled to the energy momentum tensor itself [34] . Given 
the correlation function 

C(x) = ±([* ij (x),w ij (0)]), (1) 

where irij = Tjj — is the trace-less part of the spatial components of the energy momen- 

tum tensor, the shear viscosity is the leading coefficient of its Fourier transform C(w,k) for small 
frequency and zero spatial momentum, 

n = lim — '- . 2 

First we assume that we are at zero temperature and find a simple representation for the energy- 
momentum correlation function. We will work in Euclidean space and start with the definition, 

M^ ea {x) = ( T {T^u (x)Tg a (0)) ) , (3) 

from which the following symmetry properties are evident, 

M^ v = M vfx = ag = M Qa ^ . (4) 

For the shear viscosity we will compute the following Euclidean expression: 

M W = ^£< T (^K(o))) • ( 5 ) 

ij 

At zero temperature we have 50(4) rotational invariance in Euclidean space which can be used 
to relate the value of M at a generic x^ argument to the value of Qa at r = (|x|, 0, 0, 0). 

2.1 Symmetries 

To this end M^ v e<y is decomposed into invariant tensors compatible with the symmetry properties 
(4). The most general such tensor depending on x^ can be written as, 

= M 1 (\x\) 8^8 grT + M 2 (\x\) {5 M S va + 8 ve 8^) + M 3 (\x\) X » xJp ° ± XpxJpv + 

4{\x\)— ? z-z- + M 5 (\x\) — f- — , (6) 



3 



where the coefficient functions Mi, M 2 , M3, M4 and M5 can only depend on \x\. Using this formula 
at coordinate r = (|x|, 0, 0, 0) we obtain the following special cases 



(7) 



From the above formulae one can express the functions M li2 ,3,4,5 as linear combinations of certain 
components of M^ vp(T , 



M 00 00 (t) 


= Mi(|x| 


+ 2M 2 {\x\) + 2M 3 (\x\) + 4M 4 (|x|) + M 5 (\x 


M 00 11 (t) 


- Mi(|a;| 


+ M 3 (\x\) 


Moioi(r) 


= M 2 (\x\ 


+ M 4 (\x\) 


Mn 11 (r) 


= Mi(|a;| 


+ 2M 2 (\x\) 


Mn 22 (r) 


= Mi(|a;| 





Mi (|a; 
M 2 (|a; 
M 3 (|x 

M 4 (|ar 
M 5 (|x 



) = Mn 2 2 (t) 
) = ^Mhii(t 

) =Mooii(t)-Mh22(t) 
) = Moioi(r) - ^Mlih(t) + ^Mii2 2 (r) 
) = M 000 o(t) + Mn H (r) - 2M 00 ii(r) - 4M 010 i(r). 



) = ^M U u(t)-^M 1122 (t) 



(8) 



After contracting M pu pa appropriately in order to get Ai (x) at any generic x p argument we arrive 
at, 

M(x) = M 2 (\x\) + luM A {\x\) + lu 2 M 5 (\x\), 
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where u = x 2 /x' \. This can be rewritten using (8) as 



2 4 u 

M(x) = — u 2 M ooo(t) - —u 2 Moou(t) + —(10 - 4u)M i i(t) + 

+ 4( 15 - 10u + 2« 2 )Mn 11 (r) + \{2u - 3)Mn 22 (r) . 
30 



(10) 



In the above formula the M pvpa components on the right hand side are evaluated at the point 
r = (\x\, 0,0,0). 



2.2 Plaquette-correlators 

All of the above manipulations were done in the continuum and now we will switch to the lattice by 
writing some of the energy momentum tensor components in terms of plaquette variables P pu {x) [26], 



T -I 



(11) 
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or more explicitly, 



Too — 


a 

N 


[-P01 


-P02 


-P03 


+ P21 + P31 


+ P32 


Til = 


P 
N 


[-P10 


-P12 


-Pl3 


+ P 2 + P30 


+ P32 


T 22 = 


P 
N 


[-P20 


-P21 


— P23 


+ Pw + Pso 


+ P31 


T33 = 


P 
N 


[-P30 


-P31 


— P32 


+ P20 + P10 


+ P21 



(12) 



Substituting the above back into M. we find, 
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M{x)= (15 - 10« - 4u 2 ) \2P W21 (t) - 2P W32 (t) + P ww (t) + P 2121 (t) - 

-Pi02o(r) - P 2 i3i(r)^ + ^u(5u - 2)M 0101 {r). (13) 

where P pv ga {x) = (P fll ,(x)P ea (0)) is the plaquette-plaquette correlator. 

As we will see, at small f3 i.e. large coupling, the leading and equal contributions come from 
P1021 and i-2i2i- Hence we can write 

M(x) = |^(15- 10u-4u 2 )P 2121 (t), (14) 

at leading order. 

For the shear viscosity we need the Fourier transform of M at momentum p — (po,0, 0, 0) where 
Pa -> 0, 

M W = ljf* /^(15-10«-4 U 2 ) e ^ t P 2121 (r). (15) 

The P2121 term only depends on |x| and the terms containing the u variable only depend on the 
angles. For small p one may substitute the 15 — lQu — 4u 2 term with its integral over the angles, 

M( Po ,P = 0) = |! J^xe^P^r). (16) 
Thus the plaquette-plaquette correlator is needed for plaquettes in the 21 plane. 



2.3 Strong coupling expansion 

The plaquette-plaquette correlator, G pp (t), will be calculated via the strong coupling expansion. 
The Wilson lattice formulation of pure SU(N) gauge theory is used and the action is, 

S = ^S p =/?£(l-lReTrt/ p y (17) 

P p 

where (i — 2N/g^, g is the bare coupling constant, the sum is over all plaquettes p of the lattice 
and Up is the plaquette variable built out of the links U^(x) in the usual way. 
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The strong coupling expansion amounts to an expansion in small (5; for details see [33] or 
[35]. Using the characters \R OI SU(N) there is an expansion which sums over all irreducible 
representations R, 



R=iO 



(18) 



where d R is the dimension of the representation and c(/3) and a R (P) are some coefficients. These 
coefficients will be expanded in (3. 

At leading order for the group SU(N) we have [35], 



(19) 




N > 2 



where / refers to the fundamental representation. 

The plaquette-plaquette correlator can be computed from the partition function if the coupling 
f3 is changed temporarily to /3i and /3 2 at the two plaquettes and remains f3 elsewhere 



(P 1 (t)P 2 (0)) = N 



2 lnZ(/3 1; /3 2 ,/3) 



dPid(3 2 



(20) 



/3i=/3 2 =/3 



The logarithm of the partition function is computed by forming closed connected surfaces and 
assigning representations to the plaquettes of the surface following the usual recipe of the strong cou- 
pling expansion [35] . To leading order the fundamental representation is assigned to each plaquette 
leading to the contribution 

<P f a n f (J3), (21) 

where a/(/3) is the expansion coefficient from (18) corresponding to the fundamental representation, 
df = N is its dimension and n is the number of plaquettes on the surface. The leading order terms 
in the expansion of a/(/3) were given in (19). 

The smallest surfaces spanning between well separated plaquettes are tube-like objects: they 
are closed at both end, and go almost straight from one plaquette to the other. An extra factor 
comes from the fact that the orientation of the tube and the orientation of the plaquette is not 
necessarily the same. In table 1 we summarize these extra factors for the orientations that appear 
in (13). To estimate the value of M ioi we recall that in terms of plaquettes its expression is [26], 



-P1021 


-P1032 


-P10 10 


P21 21 


Pia 20 


-P21 31 


-M0101 


1 




a f 


1 


a f 


a f 


< aj 



Table 1 : Extra factors necessary for closing the surface in case of different orientations of the ending 
plaquettes. The tube is in the temporal direction. 



Moioi(z) = (TrPoAx)Pia(x)) (Tr P 0g (0)P le (0)) ) . 



(22) 



G 



where 

Pfivix) = — - (P^ u (x) — Pnv(x)) tlaceleaspllIt ■ (23) 

The minimal surface which contains all the plaquettes of the above formula contains (in addition 
to the tube contribution) extra 2 plaquettes. Therefore the 4-plaquette expectation value has an 
extra contribution of a 2 which explains the last column of table 1. Hence we are indeed left with 
the single contribution of P21 21 as promised. 

The logarithm of the partition function in (20) at leading order is then 

lnZ(/3i,ft,/3) = d 2 f a f ((3 1 )a f ((3 2 )G E (T) , (24) 

where G e (t) is the open tube contribution, leading to 

P21,2i{t)=A n G e {t), (25) 

where 

A -N^ {P) -l l for TV = 2 

A » ~ N -p- - 1 1/4 {0lN>2 (26) 
Therefore the quantity (16) we are interested in is simply, 

M( Po ,P = 0) = ^^G E (po,p = 0). (27) 

2.4 Straight line contribution 

Now we proceed with the calculation of the open tube correlation function G e {t). 

The lowest order contribution from the surface which includes plaquette p on the hypersurface 
at t = and at t = r, both having the same spatial position, is the straight tube giving, 

G E r (r)=e-°\ (28) 

where er = — Mna/(/3) is the glueball mass. The subscript "str" stands for straight. For generic 
end points we have G E r (x) — exp(— a\x\). 
We will need its Fourier transform, 

^W = l/A e -^M = () , 2/ / +1)5/2 , (29) 

where 

12tt 2 , , 

To evaluate the shear viscosity we have to calculate the imaginary part of the retarded propa- 
gator of the corresponding energy-momentum tensor components, and up to now we calculated the 
Euclidean propagator. But the retarded and the Euclidean propagators are related as 

G E (-iuj + e,k) =«G rot (w,k), (31) 
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whereas from the retarded Green function one obtains the spectral function as 

e(*0,k) =DisciG rot Oo,k) = -2lmiG rct (k + ie,k), (32) 

feo 

where Disc stands for discontinuity. Thus the retarded Green function from (29) and (31) is 

l G ltl( k ) = (1 _ fc2 / CT 2)5/2 ' ( 33 ) 

with Landau prescription, k — > k + ie. The spectral function then reads, 

g str (k) = sgn(fc )e(fc 2 - ^ 2 ) (fc2/ J A 1)5/2 • (34) 
The spectral function exhibits a mass gap; it is zero below a, the glueball mass. 
2.5 Zig-zagging lines 

Now that we have the contribution of a straight line (tube) we will also add the contribution of 
zig-zagging lines with arbitrary number of joints. 




Figure 1: Diagrams consisting of an arbitrary number of joints. Each line is a thin tube. 
A zig-zagging line with one joint has the following contribution, 

^G^(z)G^(x - z) = [d 4 zGf tI (z)Gf tr (x-z), (35) 

a J 

which turns into G^ r (fc) 2 in Fourier space. Similarly, the contribution of zig-zagging lines with n 
joints is Gf tI (k) n+1 . Summing them all leads to, 

Gf tI (k) + Gf tI (kf + Gf tI (k) 3 + ■■■= „ - • (36) 

G%r(k) - 1 

The calculation is completely analogous to the usual self-energy summation in ordinary perturbation 
theory. The spectral function corresponding to this resummation is 

eo(fc) -sgn(fc )6(fc 2 - ( T 2 )g(fc), (37) 

where we have introduced 

2A(/c 2 /cr 2 - If/ 2 
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k 2 /c 2 

Figure 2: The spectral function at A = 0.1 (solid line) and unresummed (dashed line). The two 
curves differ only in the threshold region. The resummed spectral function exhibits a peak and a 
smooth threshold behaviour. 



The spectral function is plotted on Fig. 2.5. 

By resumming the contribution of all zig-zagging lines with an arbitrary number of joints and 
arbitrary location for these joints we have essentially resummed the contribution of all curves 
since any curve can be approximated by zig-zagging lines. This resummation changed the spectral 
function in an essential way only near the threshold: it goes to zero at the threshold at k 2 = a 2 
and exhibits a peak at 

^ma X = ^(l + A 2/5 ), JW=L (39) 

Note that it smeared out the divergence at the threshold of the straight tube case. Nonetheless 
the zig-zagging resummation really affects only the vicinity of the threshold; for example, at k 2 = 
c 2 (l + 2 A 2 / 5 ) the zig-zagging propagator is smaller than the straight line propagator only by 3%. 

The spectral function above may be interpreted as a quasi particle (with mass M = a\/\ + A 2 / 5 , 
cf. (39), and width proportional to 0A 1 / 5 ). We must, however, not forget about the presence of the 
continuum, which may alter some of the consequences coming from the quasi particle assumption. 

2.6 1-loop contribution and finite temperature 

So far we were at zero temperature where according to (38) the contribution to the viscosity is zero 
because of the presence of the mass gap. At finite temperature the finite imaginary time direction 
and periodic boundary conditions has to be taken into account leading to, 

oo 

G§(r,x)= G*(T + n/T,x). (40) 

TL— — OO 
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Figure 3: The one loop diagram contribution to the viscosity. Each line is a thin tube. 

This, however, does not change the spectral function and only introduces the usual Bose factor 
n(uj) = (exp(w/T) — into the spectral representation of the Green functions. The zero tem- 
perature spectral representation, 

oo 

iGZ(T)=aJ^e-«WQo(u), (41) 
o 

turns into 

oo 00 

]T Gq(t + n/T) = a / g («MVT-M) + e u, W ) n{u))eo{u) , (42 ) 

n= — oo q 

at finite temperature where on the right hand side the zero temperature spectral function appears. 
As a consequence the retarded Green function will remain the same as at zero temperature, G^' = 

So far we have included all single tube contributions, even the ones wrapping around the imag- 
inary time direction, but did not find any nonzero contribution below the mass gap a. This is 
expected also when we examine local modifications of the propagator, since these only modify the 
value of the mass gap [33,37]. As a hint for what types of diagrams could be relevant, we recall that 
at weak coupling the phenomenon of Landau damping is similar; it yields nonzero contribution to 
the imaginary part of the propagator below the light cone at finite temperature, while this region 
is excluded at zero temperature. The most important diagram resulting in Landau damping is a 
loop diagram, where one propagator is wrapping around the imaginary axis, while the other is not. 

We also include the effect of a 1-loop strong coupling diagram, see figure 3. In complete analogy 
to ordinary perturbation theory a self-energy can be introduced and the retarded self-energy can 
be expressed from the zero temperature one as [36] 

e-w = ^ / wf G + n{ko) ) go{k)Grip - k) • (43) 

since each of the two joints of the one-loop diagram gives a factor of a/(/3) _1 . 

Using this retarded self-energy the full one-loop corrected spectral function £>i_i 00 p(p) corre- 
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sponding to the plaqucttc-plaquette correlator is 

er-ioopW = Di f zG- loop ( P ) = Disc _ — 



Disc P0 ((iG )- 1 (p)-E ret (p)) 



(Re (iGo) _1 (p) - Rc£ rot (p)) 2 + (Im^Go)- 1 ^) - Im£ rct (p)) 2 ' 

For the viscosity the above expression is needed at zero spatial momentum and small frequency 
where Gq 1 has no imaginary part, only a non-vanishing real part. On the other hand 57 et (p) can 
be neglected to leading order in the denominator. Therefore near zero frequency we have, 

Disc T rct (n) 

ei-ioopbo « 0,p = 0) « (Rc(t g o) _ 1( o ) J )2 = (*G (0)) 2 Disc£ ret (p). (45) 

Since iG£ e *(fc = 0) = G§(k = 0) we obtain, 

ei-ioo P (po « 0, p = 0) = A 2 Disc E rot (p). (46) 

Pa 

From (43) we find 

N 4 n 4 f d 4 k { 1 \ 

^ Eretw= ^/(^ (47) 

where we also used (19) to substitute the leading order value for a/(/3). 

Using (37) at p = the following four combinations of step functions will appear in (47), 

Q(k - w k )6(po - k - Wk) , 
8(-fc - w k )6(-p + fco - Wk) , 

-e(fc -wk)e(-p + fco-w k ), (48) 

-6(-fc - ^ k )6(po - &o - Wk) ■ 

The first two is nonzero only when p > 2w k or p < — 2w k , respectively. Therefore they do not 
give contribution at small frequency. The remaining two terms yield for < po < a, 

N 4 n 4 f d 4 k 

Disc E rct (p , P = 0) = - — — / -— (n(fco) - n(k - Po)) e(fc - p - Wk)^)^ - p), (49) 



Clearly for vanishing momentum Disc Po E(0) = as it should be. The linear term in p is then, 
Disc po E ret (Po) N 4 a 4 f d 4 k ( dn(k ) \ 2 



A further simplification is possible since it is consistent with all previous approximations to 
truncate the expansion 

dn(k ) _ 1 eWT _ ^ t T 



dk T (e k «/ T - l) 2 T 



e -Ko/i +o(e- A * ' 1 ). (51) 
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at the first term because the neglected terms contribute as much as higher loop diagrams. This 
leads to, 



which in turn yields the following result for the spectral function corresponding to the plaquette- 
plaquette correlator using (46), 



2i-ioop(Po«0,p = 0) A 2 7V 4 a 4 



OO oo 

■ J dkk 2 j dk e- ko/T g 2 (k). (53) 



UJ k 

Ultimately we can express rj collecting all relations from (1), (2), (5) and (27) 

oo oo 

A 2 7V 2 a 4 



V 



2tt 3 T 



/ dkk 2 J dk a e- ka ' T f{k). (54) 



3 The entropy 

The entropy density is the derivative of the free energy density with respect to temperature. 

The free energy in the strong coupling expansion is coming from closed and connected surfaces 
similarly to the discussion in section (2.3). The first temperature dependent contribution to the free 
energy and hence the first non- vanishing contribution to the entropy comes from surfaces wrapping 
around the imaginary time direction. 

The first contribution is again from a thin straight tube line wrapping around once. The value 
of the dimensionless free energy density corresponding to such a straight tube has been calculated 
in [37] 

f = -^={aT)a s {(3) iN \ (55) 

where l/^/Az = 1 and l/\A4jv>2 = 2 comes from the fact that for N > 2 the fundamental 
representation and its complex conjugate are two different representations. The factor of 3 comes 
from the 3 possible orientation of the tube (12, 13 and 23). 

In exactly the same way as was done for the calculation of the viscosity one needs to take zig- 
zagging as well as 1-loop diagrams into account. To be able to collect these effects the free energy 
density is written as 

/ = —£j=(«T) (Gf (0) - GfL (0)) . (56) 



/An 

Now we can use the spectral representation of the propagator and write 

/ = -^Ta«/ H Q + n (fc0 ) - \ sgn( fc0 )) ,(*) . (57) 

The approximation n(k a > 0) w e~ ka l T is valid just as before since the higher loop diagrams are 
not taken into account. Within this approximation one arrives at, 

>=-^ T "'jwf°-* MTm ' (58) 
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and using the specific form of the spectral function from (37) we find 

3T« 5 



oo oo 



/ = - 2 J 3 ^_ J dkk 2 J dk e- k «/ T g(k) . (59) 



io k 

The entropy density from here is then, 

oo oo oo oo 



I dkk 2 J dk a (l + ^ e - fc «/ T g(fc) « 2 ^ 3 ^J_ | dkk 2 J dk k e~ k ^ T g{k) . (60) 

w k 



Note that the resulting entropy is independent of N 2 [37,38]. This result will be used in the next 
section to calculate the shear viscosity to entropy density ratio. 

4 The shear viscosity to entropy density ratio 

Our final goal is the ratio -q/s and using (54) and (60) from the previous sections one obtains, 

V N 2 A 2 VA^ (f(k)) T 



s 3 ( ak g(k) ) T ' 

where the following short hand notation is introduced 

d 4 k 



(61) 



{F)t = J J^®( k oMk 2 -* 2 )e- ko/T F(k). (62) 

It is worth emphasizing that all dependence on N have dropped out completely (even for finite N) 
as expected from naive counting of degrees of freedom. Both the entropy density and the shear 
viscosity are expected to scale as iV 2 and hence the ratio is expected to be independent of N in 
accordance with our formula (61). 

In order to evaluate (61) expressions of the type ( (ako) n g m (k) ) T are needed. The computation 
of these terms is greatly simplified by the fact that in the strong coupling expansion Tja is small. 

A convenient coordinate system to use is 

fc = fccosh?7, k x — fc sinh r/ sin # cos 0, k y — k sinh r\ sin 9 sin <fi, k z = k sinh r] cos 9 , 

(63) 

and the relevant quantities then are given by 

oo oo 

{(ak ) n g m {k)) T = J dkk 3 (ak) n g m {k) J d V sinh 2 r] cosh" r\ e - kcoshr >/T ( 64 ) 

a 

A lengthy but straightforward manipulation of the integral in the approximation that Tja is small 
leads to 



r— rp5/2 3/2 r / o ? 5/2 

((ako) n r(k)) T = ^^^H^" 1 J dze- z ' w ^ 



(65) 
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where w = 2 A 2 / 5 T/cr and (38) has been used for the spectral function g(k). Note that all 
n, m- independent factors in (65) will drop out from the ratio rj/s, which finally can be written as 



2T 



1 







2=*V^A»/V(^), where fM-^^ — 



2z 5 / 2 
1 + z 5 



2z 5 / 2 
1 + z 5 



(66) 



This is the first main result of our paper. The plot of f(w) - which is the rj/s ratio in the A = 1 
limit - can be seen on Fig. 4. 



0.32 



0.28 



0.26 




0.24 



0.22 



Figure 4: The rj/s ratio as a function of the rescaled temperature w for A = 1 and N = 3. The 
horizontal line represents the limiting value 0.252 which is taken as the T — > T c estimate. 

We can further analyze the expression f(w) in the small and large w regions. In the large w 
limit we can omit the exponential, and we obtain the largest contribution from the peak of the 
spectral function. As a result we obtain a w-independent form. In the small w regime, on the other 
hand, the exponential will restrict the integral to small z, i.e. the threshold region. As a result the 
ratio will behave as a power law. All together, 



5 = 7V 2 A 9 / 4 v^ < 



512 



12 5 / 4 tt VA 2 / 5 a 



IT 



2^5 



, 15 V IOtta/3 



5/2 

if T < A 2 /V 
0.056 if T > A 2 / 5 ct 



(67) 



Here high temperature of course is still limited to the confinemened phase since we are using the 
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strong coupling expansion. In particular we are not expected to reach the deconfined phase and so 
the expected increase of rj/s for increasing T above T c can not be seen here. 

Let us note that the result (67) scales as N 2 in accordance with the expectation that for T < T c 
the entropy density scales as iV° and the shear viscosity as N 2 . Also, in the high temperature limit 
the T dependence drops out. It may be in connection to the fact that at high temperature the 
relevant part of the spectrum is the peak region, which is pretty similar to a real quasi particle 
peak. The physics of the low temperature asymptotics is rather different - there the relevant part 
is the threshold regime with continuous number of excitations. 

The most robust part of our prediction is the functional dependence of the 77/s ratio on the 
temperature. Unfortunately both in the absolute magnitude of the ratio and in the temperature 
scale, the lattice spacing dependent A quantity is still present. It is, in fact, a common feature 
of the strong coupling expansion results, as it was discussed in the introduction. For the physical 
result we may use the form 2 = Kf(T/To) where K and To are the "renormalized" value of the 
prcfactor and temperature scale in (66). 

To assess K we may argue that at small A (large lattice spacing) we are safely in the applicability 
range of the strong coupling expansion, but far from the continuum limit. On the other hand at 
large A (small lattice spacing) we are near the continuum limit, but far from the applicability range 
of the strong coupling expansion. As a best compromise we may take the prefactor at A w 1, which 
suggests K = 1. 

It is even harder to assess the value of To, which is approximately the "turning point" of the 77/s 
ratio towards the low temperature asymptotics. We may argue that it is the highest temperature 
end of our calculation where we are near the continuum limit and the strong coupling expansion 
applicability range. It is then natural to expect that the To scale is determined at Ass 1 and a is 
the glueball mass near the phase transition. This consideration leads to, 

1 = { U/2 for N > 2 } X /(Wiuebai.^)). (68) 

As a consequence, from (67) and (68), we estimate the shear viscosity to entropy density ratio at 
T c to be 

V _ f 0.22 for N = 2 

s \ 0.028A 2 for N > 2 • 1 > 

This is the second main result of our paper. 

5 Discussion 

In this paper we discussed the shear viscosity to entropy density ratio in pure SU (N) gauge theories. 
The strong coupling expansion of the Euclidean lattice formulation of the theory was used. This 
expansion is known to be convergent and it yields physical results for a number of low energy 
observables [33], on the other hand its renormalization is not a fully clarified topic. This is reflected 
in our result too, since both in the magnitude of the 77/s ratio and also in the temperature scale 
a lattice spacing dependent factor was present. For a physical interpretation we used the robust 
(i.e. lattice spacing independent) prediction of the calculation, which was the functional form of 
the dependence of 77/s on the temperature and assessed the scale factors at "half way" between 
the strict strong coupling applicability range and the continuum limit. As a result we estimate the 
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r)/s ratio for gauge group SU(N) at T w T c as 0.22 for TV = 2 and 0.0287V 2 for N > 2 which is 
consistent with the conjectured lower bound 1/A-k inspired by AdS/CFT ideas. 

The temperature dependence of r]/s consists of two well separated regimes. At higher temper- 
atures (but still in the low temperature phase of the gauge theory) we see a very flat function, 
slightly rising towards smaller temperatures. This is the quasi particle regime, where we are in 
good agreement with the AdS/CFT prediction. At low temperature we arrive at a different regime, 
where the near-threshold part of the spectral function dominates, and so it is rather far from a quasi 
particle picture. The "continuum number" of peaks, forming a continuum, allow to compare to the 
constructions given in [12, 14] . In particular we observe a vanishing rj/s ratio near zero temperature. 
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